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ABSTRACT 

Lattice QCD predicts a phase transition between hadronic matter and a system of 
deconfined quarks and gluons (the Quark Gluon Plasma) at high energy densities. 
Recent results from the Brookhaven Relativistic Heavy Ion Collider (RHIC) dedi- 
cated to the study of QCD at extreme densities will be discussed and compared to 
measurements obtained at the CERN Super Proton Synchrotron (SPS). 
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1 Introduction 



Quantum Chromo Dynamics (QCD) provides, as part of the standard model, a 
very successful description of strong interaction processes involving large momentum 
transfer. However, from first principles several important aspects of QCD are still 
poorly understood. Examples are color confinement, chiral symmetry restoration 
and the structure of the vacuum. Better understanding of these concepts can be 
obtained if we are able to study quarks and gluons in a deconfined state, the so-called 
Quark Gluon Plasma (QGP). 

Such a deconfined state might be created in the laboratory in heavy-ion 
collisions at the highest energies. Theoretical guidance for this comes from Lattice 
QCD calculations. Lattice QCD predicts that at an energy density e ~ 1 GeV/fm^, 
corresponding to a temperature of about 170 MeV, the system undergoes a phase 
transition from nuclear matter to a deconfined system of quarks and gluons. 
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Figure 1: a) Energy density divided by versus T at fiB = b) Transition 
temperature as a function of fiB i^l The dotted line illustrates the rapid crossover 
while the solid line illustrates the first order phase transition. 



Figure ^ shows the energy density divided by the fourth power of the 
temperature, T, versus the temperature from a lattice calculation pp. This figure 
shows that in between 150-200 MeV the energy density increases rapidly which is 
indicative of a phase transition where at high temperature the quarks and gluons 
become the relevant degrees of freedom. The figure also indicates were according 
to our current understanding the different heavy-ion machines are located on this 
diagram. These calculations are done with zero baryon chemical potential, fiB, re- 
flecting the conditions of the early universe. Small values of the chemical potential 
are obtained at RHIC collider energies whereas at lower energies, e.g. AGS and 
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SPS, the value of /i^ is large. In Fig. ^ the relation between the transition temper- 
ature and the chemical potential from recent lattice calculations 2| is shown. The 
calculation indicates that the transition temperature decreases with increasing fiB 
and furthermore that at low /i^ the transition from the hadronic phase to the QGP 
is a rapid crossover (dotted line) while at large fiB a first order transition should 
take place (full line). 

2 The Relativistic Heavy Ion Collider (RHIC) 

Heavy-ion physics entered a new era with the advent of the Relativistic Heavy Ion 
Collider (RHIC) at Brookhaven National Laboratory. RHIC is a versatile collider 
providing collisions with different ion species (ranging from protons to gold) at a 
wide range of center of mass energies ^/s^. In the three years of operation collisions 
were provided for Au+Au at 19.7, 130 and 200 GeV, p+p at 200 GeV and d+Au 
at 200 GeV. Note that the top center of mass energy for p+p is 500 GeV at RHIC. 
For details on the detectors (BRAHMS, PHENIX, PHOBOS, STAR), see 0. 

3 Event Characterization 

Heavy ions are extended objects and the system created in a head-on collision is 
different from that in a peripheral collision. Therefore, collisions are categorized 
by their centrality. Theoretically the centrality is characterized by the impact pa- 
rameter b which, however, is not a direct observable. Experimentally, the collision 
centrality can be inferred from the measured particle multiplicities if one assumes 
that this multiplicity is a monotonic function of b. Another way to determine the 
event centrality is to measure the energy carried by the spectator nucleons (which 
do not participate in the reaction) with Zero Degree Calorimetry (ZDC). A large 
(small) signal in the ZDCs thus indicates a peripheral (central) collision. 

Two other measures of the centrality which are often used are the number 
of wounded nucleons and the equivalent number of binary collisions. These numbers 
are related to the impact parameter b using a realistic description of the nuclear 
geometry in a Glauber calculation, see Fig. |21 Phenomenologically it is found that 
soft particle production scales with the number of participating nucleons whereas 
hard processes scale with the number of binary collisions. 
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4 how-pt Observables 



Examples of global observables which provide important information about the cre- 
ated system are the particle multiplicity and the transverse energy. Figure El shows 
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Figure 2: Number of wounded nucleons 
and binary collisions versus impact pa- 
rameter. 
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Figure 3: Transverse energy as a 
function of centrality as measured by 
PHENIX g]/. 



the transverse energy versus the collision centrality as measured at ^/s^ = 130 GeV 
by the PHENIX collaboration This measurement allows for an estimate of the 
energy density as proposed by Bjorken [S] 

1 1 dEr 
T^B? ctq dy ' 

were R is the nuclear radius and tq is the effective thermalization time (0.2-1.0 fm/c). 
From the measured {dET/drj) = 503 ± 2 GeV it follows that e is about 5 GeV/fm^ 
at RHIC. This is much larger than the critical energy density of 1 GeV/fm^ from 
Lattice QCD (see Fig. ID). 

Figure 0] shows the charged particle multiplicity distributions versus the 
pseudorapidity rj measured by PHOBOS at three different energies 0. The gross 
features of the particle multiplicity distributions are described by a similar behavior 
of the tails (limiting fragmentation) and a plateau at mid-rapidity consistent with a 
boost invariant region of Ay ^ 1. Notice that in total about 5000 charged particles 
are produced in the most central Au+Au collisions at the top RHIC energy. 



4 





Z 

■n 

600 
400 
200 



1 , 1 , 1 1 1 1 1 1 , 

"^"^ 19.6 GeV 
<= 6-1 5% 

.15-25% 

L □ 25-35% J 
f 35-45% 
45-55% 


1 , , 1 , 1 1 1 1 
L 130 GeV J 

^ _ 
fl • 

• ^oooooooooooooooooo, *, 

- 


1 1 1 , , 1 , 1 1 
L 200 GeV J 

/ 

• 0° °o ■ 

, - 

- '/ V - 




•.0° 

. ■■ ■■■ ° • 
■ h • 


. 0. 

,0 .■ ■■ 0, 




^^^^^^^^^^^^^^^^^^^^^^^ 


- .0°..>*^=™^".V _ 




. . , 1 . , . 1"^^ 



Figure 4: Multiplicity versus pseudo-rapidity for 19.6, 130 and 200 GeV measured 
by PHOBOS 10. 



4.1 Particle Yields 

The integrated yield of the various particle species provides information on the 
production mechanism and the subsequent inelastic collisions. A very successful 
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Figure 5: Particle yield ratios at RHIC compared with a thermal model J^. 



description of the relative particle yields is given by the thermal model. In Fig. El 
the particle yield ratios measured at RHIC are plotted and compared to values 
from a thermal model fit [H]. The results from the fit show that all particles ratios 
are consistent with a single temperature and single chemical potential in a thermal 
description. The temperature obtained in this way, 176 MeV, is called the chemical 
freeze-out temperature and the measured value is very close to the phase boundary 
value in lattice calculations (see Fig. ^p). Figure shows the relative particle 
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Figure 6: a) Particle ratios versus rapidity measured by BRAHMS b) The ratio 
K^/K^ versus p/p or, equivalently, fj,B- 

ratios of pions, kaons and protons and their anti-particles versus rapidity [Zj. For 
the heavier particles the ratio drops rapidly for y > 1. Figure IHb shows the ratio 
of K~ /K^ versus p/p for AGS to RHIC energies. The decreasing ratio of p/p as 
a function of rapidity can thus be understood from the changing baryon chemical 
potential at a constant chemical freeze-out temperature. 

4.2 Spectra 

The particle spectra provide much more information than the integrated particle 
yields alone. The particle yield as a function of transverse momentum reveal the 
dynamics of the collision, characterized by the temperature and transverse flow 
velocity of the system at kinetic freeze-out. Kinetic freeze-out corresponds to the 
final stage of the collision when the system becomes so dilute that all interactions 
between the particles cease to exist so that the momentum distributions do not 
change anymore. Figure [7^ and b show the transverse momentum distributions at 
^/s = 17 GeV from NA49 [H]. The lines are a fit to the particle spectra with a 
hydrodynamically inspired model (blast wave). The fit describes all the particle 
spectra rather well which shows that these spectra can be characterized by the two 
parameters of the model: a single kinetic freeze-out temperature and a common 
transverse flow velocity. Figure ^ shows the combined pion, kaon and proton pt- 
spectra from the four RHIC experiments. Also at these energies it follows from a 
common fit to all the spectra that the system seems to freeze-out with at a similar 
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Figure 7: NA49 (SPS) and RHIC low-pt spectra O O [HI El/. 



temperature and with a transverse flow velocity as observed at SPS energies. 
4.3 Azimuthal Correlations with the Reaction Plane 

The nuclear overlap region in non-central collisions has an almond like shape with 
its longer axis perpendicular to the reaction plane (the plane defined by the beam 
axis and the impact parameter). This particular shape leads to a pressure gradient 
which is different in and out of the reaction plane which, in turn leads to azimuthally 
asymmetric particle emission. The asymmetry can be described by: 

E—- = — - 1 + J2 2vnCos{n^)] 

d-^p 2Tiptdptdy 

where is the azimuthal angle with respect to the reaction plane and the coefficient 
of the second harmonic, f2, is called elliptic fiow. The magnitude of V2 and its pt 
dependence allows for the extraction of the kinetic freeze-out temperature and the 
transverse fiow velocity as function of emission angle. 

Figure |H1 shows the measurement of f 2 versus pt for pions and protons 
plus antiprotons. Due to transverse fiow the pt dependence of V2 depends on the 
particle mass as is evident from Fig. |H| Also shown in this figure are hydrodynamical 
model calculations using two different equations of state (THj corresponding a hadron 
gas and a QGP. It is seen that the QGP EOS shows the best agreement with 
the data. In Fig. El RHIC data on V2{pt) JHl 113 UHl HHI for various particles are 
compared to a hydrodynamical inspired blast wave fit. The agreement of the data 
with this fit shows that the V2{pt) for all particles can be described in terms of a 
single temperature and a 0-dependent transverse fiow velocity. Furthermore, the 
magnitude and pt dependence of the elliptic fiow for the various particles suggest 
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Figure 8: V2(pt) for pions and protons at 
y/s^ =130 116]. The lines are hydro- 
dynamical model calculations 11^ . 



Figure 9: V2(pt) for identified particles 
compared to a blast wave fit f73|. \F\ [7^ 



strong partonic interactions in an early stage of the collision and, perhaps, early 
thermalization of the system. 

5 High-p( observables 

In heavy-ion collisions at RHIC, jets with transverse energies above 40 GeV are pro- 
duced in abundance, providing a detailed probe of the created system. However the 
abundant soft particle production in heavy-ion collisions tends to obscure the char- 
acteristic jet structures. At sufficient high-p^ the contribution from the tails of the 
soft particle production becomes negligible and jets can be identified by their leading 
particles. It was proposed that a leading particle traversing a dense system would 
lose energy by induced gluon radiation (so called jet-quenching [201 )■ The amount 
of energy loss is in this picture directly related to the parton density (mainly gluons 
at RHIC) of the created system. Currently there are three observables sensitive to 
this energy loss as discussed in the next two subsections. 

5.1 Single Inclusive Particle Yields 

As mentioned above, the single inclusive particle yield at sufficiently high-pf is domi- 
nated by the leading particles from jets. Figure ITUk shows the ttq spectra as measured 
in p-|-p at ^/s = 200 GeV. In the same figure also two NLO QCD calculations are 
shown. The ratio of the data to the theory shows that in p+p the ttq spectrum 
is well described. In Fig. ITUfe the charged hadron spectra measured in Au+Au at 
y/s^ = 200 GeV and the p-|-p reference spectra at the same energy are shown. 
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Figure 10: ttq spectra in p+p and charged hadrons in p+p and Au+Au Jd^. 



One of the observables suggested for measuring energy loss is the so called nuclear 
modification factor defined by 



RAA^Pt 



dVAA/dydpt 



(A^bmary)dVpp/d?/dpt' 

where d'^app/dydpt is the inclusive cross section measured in p+p colhsions and 
(iVbinary) accouuts for the geometrical scaling from p+p to nuclear collisions. In 
the case that a Au+Au collision is an incoherent superposition of p+p collisions 
this ratio Raa would be unity. Energy loss and shadowing would reduce this ratio 
below unity while anti-shadowing and the Cronin effect would lead to a value above 
unity. Figure ITTb.c shows this ratio for charged particles and tto's in central Au+Au 
collisions at mid-rapidity. The ratio is well below one and at high-p^ the suppression 
is a factor of 5. At intermediate pt the charged particles and ttq are both suppressed; 
however the magnitude differs by a factor of two. In Figure ITTH Raa is plotted at 
more forward rapidities showing that the suppression also persists there. 

To discriminate between energy loss and shadowing, d+Au collisions were 
measured. If the suppression is due to shadowing it should also be observed in the 
d-Au system. Figure ITTk shows the d+Au spectra versus centrality and Fig. ITTb.c 
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Figure 11: d+AU and Au+Au measurements from PHOBOS STAR 
PHENIX IMJ and BRAHMS JE^. 



the nuclear modification factor for charged particles and ttq, respectively. It is clear 
that in d+Au interactions no suppression is observed. In fact, to the contrary, a 
small enhancement is seen consistent with the Cronin effect. From this observation 
it follows that the observed suppression in Au+Au collisions is due to final state 
interactions. The magnitude of the observed suppression at the top RHIC energy 
indicates, in the jet quenching picture, densities which are a factor 30 higher than 
in nuclear matter. 

5.2 Azimuthal Correlations 

In heavy-ion collisions, azimuthal correlations between particles can be used to study 
the effect of jet quenching in greater detail. The azimuthal correlations of two high-p^ 
particles from jets are expected to show a narrow near-side correlation and a broader 
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away-side correlation. However, in the case of strong jet quenching the away-side 
jet would suffer significant energy loss and would be suppressed. Recently, CERES 




Figure 12: Azimuthal correlations measured at the SPS by CERES \2f^ . 



measured such a correlation function at the top SPS energy. In Fig. IT^ the near- 
side correlation (at A</) = 0) shows a narrow peak consistent with the correlation 
observed in jets. The away-side correlation peak is observed in more peripheral 
collisions but disappears for more central collisions, see Fig. 112b. Figure IT^ shows 
that the width (a) of the near-side correlation peak stays constant as a function of 
centrality, but that the away-side peak broadens for more central collisions. The 
total integrated yield is the same in the near and away-side peak (Fig. ^Ji). There- 
fore, the disappearance of the away-side peak at the top SPS energy is interpreted 
as being due to initial state broadening 




Figure 13: a,h) Back to hack correlations \27l and c) elliptic flow parameter at 
intermediate pt- 

The azimuthal correlations of high-pf particles (trigger particle 4 < pt < 
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6 GeV/c, associated particle 2 GeV/c < pt < pP^) measured in p+p collisions at 
RHIC are shown as the histogram in Figure fT^. The near-side and away-side peaks 
are clearly visable. The correlation function observed in central Au+Au collisions 
(stars in Fig. {T^o) shows a similar near side peak while the away-side peak has 
disappeared. 

To investigate if this is due to initial state effects, the same analysis was 
done for d+Au collisions. In Fig. IT^ the near and away-side peaks are shown 
for minimum bias and central d+Au collisions compared to p-|-p. The away-side 
correlation in d+Au is clearly observed even for the most central collisions. Com- 
paring the away-side correlation in p+p, d+Au and Au+Au, Fig. IT^ . shows that 
the suppression only occurs in Au+Au collisions and therfore is a final state effect 
as expected from jet quenching. 

The energy loss depends on the distance traversed through the dense 
medium by the partons. In a non-central collision the distance will depend on 
the azimuthal angle with respect to the reaction plane (see \ow-pt section). Because 
the hard scattering producing the di-jet has no correlation with the reaction plane, 
an observed asymmetry in the high-pf particle emission will be due to final state 
interactions (such as the jet quenching mechanism). In Fig. 113b the observed elliptic 
flow signal as a function of pt is shown for charged particles, kaons and lamdas. It 
is clear from this figure that the observed asymmetry is very large up to the highest 
Pt measured. Like the nuclear suppression factor Raai the elliptic flow at interme- 
diate Pt depends on the particle species. This could be due to an interplay between 
the soft hydrodynamical behavior and the jet quenching, which would cause a mass 
dependence. However, more recently this has been interpreted as a possible sign 
of particle production at intermediate pt by parton coalescence. In that case it is 
not the mass of the particle which is responsible for the splitting but rather the 
number of constituent quarks (two for mesons and three for baryons). A deflnitive 
test will be the measurement of elliptic flow of the (/)-meson because in the coales- 
cence interpretation it should have an elliptic flow similar to the pions while in the 
hydrodynamical interpretation is would have an elliptic flow value similar to the 
proton. 

6 Conclusions 

The flrst three years of RHIC operation have provided a wealth of interesting data. 
We have seen that: 

• Particle yields indicate a chemical freeze-out of the system near the phase 
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boundary; 



• Identified particle spectra are consistent with boosted thermal distributions 
and identified particle elliptic fiow shows remarkable agreement with ideal 
hydrodynamical calculations based on a QGP equation of state; 

• The particle yield at high-pt is suppressed compared to proton-proton reference 
data. The fact that this suppression does not occur in d+Au collisions shows 
that it is a final state effect, consistent with parton energy loss in dense matter 
(jet quenching); 

• The suppression at intermediate pt shows a particle dependence which could 
be explained by particle production, at intermediate pt, by parton coalescence; 

• The elliptic fiow at intermediate pt is large and also shows a particle depen- 
dence. Like above, this is consistent with energy loss in dense matter and 
particle production via parton coalescence; 

• In the most central events the high-p^ back to back correlations are consistent 
with zero. Such disappearance of the away-side jets is expected in the case of 
very strong energy loss in a dense medium. 

All these observations, taken together, are consistent with the creation 
of a very dense and strongly interacting system in heavy-ion collisions at RHIC 
energies. While all these observations are consistent with the creation of a QGP, 
more detailed knowledge of QCD at high densities and temperatures is required. 
This poses a formidable challenge for theory but will be crucial for the detailed 
interpretation of the present and future data taken at RHIC and LHC. 
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